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THE  ELECTROMAGNETIC  BIAS  OF  ALTIMETER  MEASUREMENTS 
OF  MEAN  SEA  LEVEL  AS  DETERMINED  BY  AN  AIRBORNE 

10  GHz  RADAR 


INTRODUCTION 


Satellite  radar  altimetry  has  matured  greatly  since  the  S-193  microwave 
experiment  on  SKYLAB  ten  years  ago.  For  geodesy  and  oceanography,  accurate 
radar  altimeter  data  are  needed  for  mapping  of  the  marine  geoid,  ocean 
currents,  and  mesoscale  features.  While  the  radar  altimeter  is  a  range 
measuring  device,  ocean  wave  heights  and  wind  speeds  at  the  sea  surface  are 
also  obtained  with  proper  analysis  of  the  data.  However,  measurement  of 
range  to  high  accuracy  in  the  presence  of  ocean  waves  is  not  problem  free. 
The  diameter  of  the  footprint  illuminated  by  a  satellite-borne  radar 
altimeter  is  large  compared  with  the  lengths  of  the  ocean  wind  waves.  For 
example,  the  pulse-width-limited  footprint  of  the  SEASAT  altimeter  had  a 
diameter  of  1.6  km  for  smooth  seas.  The  large  footprint  effectively  acts  as 
a  spatial  filter  (Figure  l)  that  averages  the  height  of  many  waves;  thus  the 
amplitudes  of  ocean  waves  do  not  directly  enter  into  the  range  measurement. 
However,  it  has  been  observed  (Yaplee,  et  al,  1971;  Shapiro,  et  al,  1972; 
Kenney  and  Walsh,  1978;  Choy  and  Uliana,  1980)  that  wind  generated  ocean 
wave  troughs  tend  to  reflect  nadir  incident  microwave  signals  better  than 
the  wave  crests.  Since  a  radar  altimeter  determines  its  range  by  measuring 
the  round  trip  time  of  a  pulse  of  electromagnetic  energy  to  the  sea  surface, 
and  uses  the  centroid  of  the  pulse  as  the  time  reference,  the  range  so 
derived  tends  to  be  biased  toward  the  troughs.  If  we  define  mean  sea 
surface  as  the  centroid  of  the  wave  height  distribution,  we  see  that  the 
radar  determined  "mean"  sea  surface  may  not  be  the  same.  The  difference  has 
been  called  the  electromagnetic  bias  or  the  EM  bias.  To  achieve  a  range 
accuracy  of  a  few  centimeters  with  spaceborne  altimetry,  understanding  the 
EM  bias  problem  is  important. 

The  first  experiment  to  investigate  the  EM  bias  was  reported  by  Yaplee, 
et  al  (1971).  Their  radar  measurements  were  taken  from  a  navigation  tower 
off  the  coast  near  Norfolk,  Virginia.  Only  two  cases  were  reported  with  low 
to  moderate  wave  heights.  Later  reports  on  EM  bias  used  data  collected  by 
airborne  sensors  with  a  wider  range  of  wave  heights  (Kenney,  et  al,  1979; 
Choy  and  Uliana,  1980).  There  have  also  been  reports  on  estimates  of  EM 
bias  using  data  from  satellite-borne  radar  altimeters  (Lipa  and  Barrick, 

1980 ;  Born,  et  al,  1982;).  Theoretical  ocean  wave  models  which  explain  the 
EM  bias  have  also  been  reported  (Jackson,  1979;  Lipa  and  Barrick,  1980). 
Comparisons  among  various  observations  and  observations  with  theories 
indicate  that  the  EM  bias  question  is  not  completely  resolved.  The  present 
investigation  is  part  of  a  Joint  program  between  NRL  and  NASA  (GSFC/WFC)  to 
address  this  EM  bias  problem. 
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RADAR  ALTIMETER  GEOMETRY  ILLUSTRATING  THE  SPATIAL  FILTERING 


H  =  HEIGHT  OF  ALTIMETER 

ii  =  ANGULAR  BEAM  WIDTH  OF  ALTIMETER 

FP  =  ALTIMETER  FOOTPRINT  OR  SPOT  SIZE  ON  THE 
OCEAN  SURFACE 

Figure  1  -  Radar  Altimeter  Geometry  Illustrating  Spatial  Filtering 


THEORETICAL  BASIS  FOR  THE  EXPERIMENT 


* 

The  airborne  radar  altimeter  views  the  ocean  in  the  nadir  direction  and 
measures  the  ocean  surface  elevation  profile  and  the  backscattered  radar 
power.  From  these  measurements,  one  can  construct  the  probability  density 
function  of  the  ocean  surface  elevations  and  the  radar  surface  impulse 
response,  both  quantities  as  functions  of  vertical  ocean  wave  displacement 
relative  to  the  displacement  mean.  The  difference  between  the  mean  values 
of  the  probability  density  function  of  the  ocean  surface  elevations  and  the 
radar  surface  impulse  response  is  defined  as  the  electromagnetic  range  bias 
(EM  bias)  of  the  altimeter.  The  radar  energy  is  confined  to  a  narrow  pulse 
(l  ns)  illuminating  the  ocean  surface  with  a  footprint  which  is  small 
relative  to  the  dominant  ocean  wavelengths  encountered.  A  software  radar 
altimeter  range  tracker  is  constructed  to  track  on  the  centroid  or  the  mean 
value  of  the  return  power.  Thus,  the  probability  density  of  the  surface 
elevations  can  be  approximated  by  the  ratio  of  the  number  of  pulses  received 
for  a  given  range  resolution  cell  to  the  total  pulses  received,  and  the 
radar  surface  impulse  response  is  approximated  by  the  ratio  of  the  average 
backscattered  power  in  a  range  cell  to  the  total  power  received  (Figure  2). 

Higher  moments  of  the  wave  elevation  distribution  (rms  wave  height, 
skewness,  and  kurtosis)  can  be  obtained  from  the  probability  density 
function  after  carefully  removing  the  contamination  of  the  aircraft  motion 
from  the  data.  However,  due  to  the  finite  illuminated  footprint  of  the 
radar  altimeter,  the  values  of  the  moments  are  underestimated  by  the  spatial 
filtering  effect.  The  numerical  reduction  in  the  estimated  values  of  the 
variance  and  skewness  of  the  wave  height  distribution  can  be  calculated  to  a 
first  order  by  using  a  simple  Phillips  spectral  model  F(K)  for  the  water 
waves  (Jackson,  1979). 


u20 

F(K)  =  eK“3  K  <  K  <  » 

o 

where 

Kq  s  pif /L  is  the  dominant  wave  number  and  L  is  the  dominant  wavelength. 
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OCEAN  SURFACE  ELE 


(DAR  RETURN  POWER  (ARBITRARY  UNIT) 


»2q  is  the  variance  of  the  wave  height  probability  density, 
vi  2q  i-8  the  third  moment  of  the  wave  height  probability  density, 
is  the  normalized  skewness. 

nKQ  represents  the  high  frequency  cut  off  due  to  the  radar  footprint 
diameter  and  is  equal  to  2k  /D. 

D  is  the  footprint  diameter. 

Carrying  out  the  indicated  integration,  the  spatially  filtered  values  are 


It  can  be  seen  that  when  n  is  large,  more  o; 
contribute  and  both  and  u ' „  approach  8/2K 

tered  values  of  UgQ  ana  respectively.  ° 


the  high  frequency  waves 

and  28/K  ,  the  unfil- 
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The  upper  limit  of  integration  is  equated  with  the  wave  number  whose 
wavelength  equals  the  footprint  diameter  D.  It  follows  that  n  is  the  ratio 
of  dominant  wavelength  to  footprint  diameter. 


n2ir 

L 


n 


Since 


X 


30 


3/2 


and  the  filtered 


or 


Table  1  lists  the  values  of  M*  20^20  an<*  X*30^X30  ^0r  severa^-  valu®8 
of  n.  ~  ~ 


Table  1 


The  Values  of  ^20^20  an<i  ^30^30  aS  func’tion  °"  n 
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The  NRL  radar  altimeter  has  a  footprint  diameter  about  15  meters  for 
the  data  used  in  this  report.  Thus  for  ocean  surfaces  with  dominant  wave¬ 
lengths  of  75  meters  or  greater,  the  reduction  in  the  estimates  of  variance 
and  normalized  skewness  is  k  percent  or  less. 

EXPERIMENTAL  PROCEDURES 

The  high  resolution  radar  altimeter  antenna  and  the  laser  profilometer 
were  mounted  near  to  each  other  on  an  airborne  platform.  The  illumination 
spot  of  the  radar  is  larger  than  the  laser  by  a  factor  of  50  at  all  alti¬ 
tudes.  Both  sensors  were  aligned  to  look  in  the  nadir  direction  under 
stable  flying  conditions.  The  sampling  rate  of  the  laser  is  200  per  second. 
The  sampling  rate  of  the  radar  altimeter  can  be  set  on  39  or  78  samples  per 
second.  The  aircraft  ground  speed  was  approximately  100  meters  per  second. 
Thus  tbs  range  profile  produced  simultaneously  by  both  sensors  are  along  the 
same  sea  surface  track.  The  laser,  in  addition  to  providing  an  independent 
check  on  the  performance  of  the  radar,  offers  a  quantitative  look  at  the 
spatial  filtering  effect  of  the  radar.  The  NRL  radar  is  a  10  GHz^  bistatic, 
CW  pseudo-random  noise  coded  system.  The  antenna  beam  width  of  6  combined 
with  the  one  nanosecond  range  resolution  assures  pulse  width  limited 
illumination.  The  laser  profilometer  is  an  IR  pulsed  laser  with  a  pulse 
width  of  10  ns. 

The  NRL  airborne  platform  is  a  Navy  P-3  aircraft.  It  is  equipped  with 
an  inertial  navigation  system  which  includes  the  monitoring  of  wind  speed 
and  direction.  A  separate  set  of  sensors  are  installed  in  the  laser  profil¬ 
ometer  housing  close  to  the  radar  altimeter  to  record  pitch,  roll,  and 
vertical  acceleration  information  of  the  aircraft.  For  most  of  the  data 
examined,  the  aircraft  motion  changed  slowly.  Of  the  three  motions,  pitch, 
roll,  and  vertical  acceleration,  the  vertical  motion  was  dominant.  Typi¬ 
cally,  in  a  three-minute  interval,  two  or  three  large  changes  in  vertical 
velocity  (approximately  50  cm/ sec)  can  be  observed.  During  stable  flight 
segments  vertical  velocity  of  the  aircraft  was  less  than  5  cm/sec.  The 
aircraft  vertical  motion  was  removed  by  an  algorithm  based  on  double 
integration  of  the  accelerometer  data  and  high  pass  filtering. 
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All  data  used  in  this  report  were  taken  at  altitudes  of  500  and  750 
feet  in  three  locations,  the  North  Sea,  Gulf  of  Alaska,  and  the  Atlantic 
Ocean  just  off  the  East  Coast  of  the  United  States.  When  the  aircraft 
reached  the  operating  altitude,  data  were  taken  continuously  in  3  to  10 
minute  intervals  and  recorded  on  digital  magnetic  tape.  During  one  of  the 
data  flights  reported  in  this  paper  the  NRL  and  NASA/Wallops  radar  systems 
were  flown  on  the  same  aircraft.  On  another  data  flight  reported  here  the 
two  radar  systems,  although  mounted  in  separate  airplanes,  collected  data 
simultaneously  at  the  same  location.  These  experiments  conducted  with 
NASA/Wallops  were  at  locations  off  the  east  coast  of  the  United  States. 

DATA  ANALYSIS  PROCEDURE 

Data  records  were  divided  into  two-minute  segments.  Altimeter  height 
measurements  were  computed  using  a  software  range  tracker  which  operates  on 
the  centroid  of  the  averaged  (4  samples)  received  radar  pulse  waveform.  The 
backscattered  power  was  determined  for  each  pulse  using  the  total  ->ower 
contained  in  the  averaged  ( U  samples)  waveform.  Although  other  t  of 
range  trackers  using  either  the  location  of  the  peak  amplitude  c  ;he 

waveform  or  a  leading  edge  threshold  value  on  the  waveform  were  t  1 ,  we 
settled  on  the  centroid  tracker  because  the  result  obtained  with  was  more 
stable.  The  first  two  moments  of  the  wave  height  distribution  we  'imputed 
from  height  data  using  standard  methods. 

N  N  2 

woo  =  TT.^  hi  u20  =  N^l”.^^hi  "  uoo^ 

i=l  i=l 


u  =  first  moment 
Hoo 

h^  =  sample  of  height 
N  =  number  of  sample  heights 


UgQ  =  second  moment 

The  third  and  fourth  moments  (skewness  and  kurtosis)  were  estimated  using  an 
ocean  surface  elevation  density  model,  Longuet-Higgins  ( 19^3 ) • 

_  _  1  /  T^\l  l\  .  X  Lz  -,m\  .  K  Lh  .  AT 


where 


-k  1  *  ’j 


g  =  probability  density  of  height. 


T  =  t/<rs. 

t  =  the  height  above  the  mean  sea  level. 

<r  =  the  ms  wave  height, 
s 

A  =  the  skewness  of  the  wave  height  distribution. 
K  =  the  kurtosis  of  the  distribution. 
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The  model  is  a  low  order  case  of  a  general  probability  function  used  by 
Longuet -Higgins  (19^3)  for  a  random  variable  that  is  weakly  nonlinear. 

For  a  given  a  ,  parameters  \  and  K  were  varied  so  the  model  best  fits 
the  data  in  a  leas?  square  sense.  The  conventional  Chi-square  is  used  to 
estimate  the  goodness  of  fit.  The  computation  of  skewness  is  extremely 
sensitive  to  the  data  noise  in  the  tail  ends  of  the  wave  height  distribu¬ 
tion.  Data  which  gave  a  poor  fit  to  the  model  also  gave  unreasonable 
estimates  of  skewness  if  the  usual  formula  from  statistics  was  used. 

For  the  majority  of  the  observations,  the  radar  sampling  rate  was  set 
at  39 /second.  For  a  two-minute  segment,  the  number  of  data  points  that  went 
into  each  determination  of  EM  bias,  significant  wave  height,  and  skewness  of 
the  wave  height  distribution  is  over  4000.  The  dominant  wavelength  informa¬ 
tion  was  derived  by  taking  a  one-dimensional  power  spectrum  of  the  wave 
profiling  data,  flying  in  the  upwind  or  downwind  direction.  Since  we  have 
no  other  independent  method  of  determining  the  wave  propogation  direction, 
the  dominant  wavelength  values  were  rough  estimates.  The  wind  speed  was 
recorded  at  the  aircraft  altitude  of  500  ft  and  no  attempt  was  made  to 
re-scale  the  wind  speed  to  a  value  corresponding  to  the  sea  surface. 

DISCUSSION  AND  CONCLUSIONS 

For  significant  wave  heights  (SWH)  up  to  five  meters,  the  value  of  EM 
bias  increases  in  proportion  to  SWH.  To  facilitate  the  comparison  of  EM 
bias  with  other  ocean  surface  conditions,  the  unit  of  EM  bias  will  be 
expressed  as  a  percentage  of  the  SWH  for  this  section  of  the  report. 

Figure  3  is  a  plot  of  EM  bias  (percentage  SWH)  mr s.  SWH.  Data  that  were 
taken  on  the  same  day  are  represented  by  the  same  symbol.  For  most  of  the 
data  examined  with  significant  wave  height  in  1  to  5  m  range,  the  EM  bias  is 
between  three  to  five  percent  of  SWH.  On  one  occasion,  however,  the  EM  bias 
was  centered  at  1.5%  of  the  SWH.  On  that  day,  the  NRL  radar  was  flying 
jointly  with  the  NASA/Wallops  experiment  on  the  NASA  aircraft.  The  skewness 
of  the  wave  height  distribution  was  observed  to  be  rather  small  indicating 
that  the  sea  was  dominated  by  swell.  Since  the  data  are  grouped  by  the  day 
of  observation,  the  error  bars  were  placed  only  on  one  representative  member 
of  a  data  group.  In  general,  the  EM  bias  estimates  have  a  standard  devia¬ 
tion  of  about  +  0.5  percent  of  SWH.  The  SWH  estimates  have  standard 
deviations  of  about  +  0.15  x  SWH  in  meters.  These  uncertainties  are  due  to 
the  limited  number  of  data  samples  available. 

Figure  4  is  a  plot  of  the  EM  bias  vs.  the  skewness  of  the  wave  height 
distribution.  It  is  of  interest  to  compare  the  regression  line  of  these 
data  with  the  result  of  Jackson  (1979)*  Jackson  concluded  that  the  EM  bias 
expressed  in  meters  is  the  product  of  the  skewness  times  the  rms  wave 
height.  This  relation  is  shown  by  the  straight  line  in  the  figure.  The 
line  has  a  slope  of  0.04  and  an  intercept  of  zero.  The  regression  line  for 
the  data  has  a  slope  of  0.03  and  an  intercept  of  0.03.  Again  only  a  few 
representative  error  bars  are  shown.  These  errors  were  calculated  using  the 
relation 
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r  is  the  mean  squared  a  moment  of  the  estimation  error  using  N 

3* 

samples , 

u  is  the  a  moment  of  the  data  samples, 
v  is  the  2a^  moment  of  the  data  samples,  and 
N  is  the  number  of  independent  samples. 

The  number  of  data  samples  were  on  the  order  of  4000  for  every  point. 
The  errors  on  the  estimates  can  be  obtained  by  assuming  that  the 
distribution  of  the  surface  elevations  is  Gaussian  with  zero  mean  and  a 
standard  deviation  as  obtained  from  the  profiles  of  the  surface.  Figure  2 
shows  a  typical  wave  height  distribution  and  backscattered  power  distribu¬ 
tion  as  observed  by  the  radar. 

To  test  for  any  dependance  of  the  EM  bias  (%  SWH)  on  several  surface 
conditions,  all  of  the  63  combinations  of  SWH,  skewness,  kurtosis  dominant 
wavelength,  wind  speed,  and  significant  slope  (Huang  and  Long,  1980)  were 
used  to  fit  the  measurements  to  multidimensional  straight  lines.  For  each 
combination  the  standard  deviation  of  the  measured  values  of  EM  bias  (%  SWH) 
relative  to  the  best  fit  line  was  computed.  Figure  5  shows  the  standard 
deviations  of  the  difference  between  the  measured  values  of  EM  bias  and  the 
corresponding  least  squares  value  displayed  in  the  order  of  the  largest 
value  to  the  smallest  value.  The  poorest  fit  resulted  when  using  SWH  alone 
(S.D.  *  1.2%  SWH)  and  the  best  fit  was  found  when  all  conditions  were 
included  (S.D.  *  0.5%  SWH).  The  dots  along  the  bottom  of  Figure  5  show  by 
their  position  the  surface  conditions  used  for  a  particular  combination. 

For  example  in  Figure  5  the  line  that  joins  the  two  dots  in  the  rows  of  wind 
speed  and  significant  slope,  showing  that  these  two  conditions  were  used, 
are  located  directly  below  an  upper  dot  indicating  that  the  standard  devia¬ 
tion  of  the  fit  was  0.65%  SWH.  There  doesn't  appear  to  be  any  pronounced 
improvement  in  the  fit  for  any  combination  of  the  conditions. 

Future  satellite  altimeter  programs  need  to  consider  the  EM  bias  as  a 
source  of  error.  Care  should  be  exercised  in  assigning  any  particular 
value  to  the  EM  bias.  Our  data  shows  a  trend  in  the  direction  of  Jackson's 
analytical  result.  Our  result  together  with  (Walsh,  1982  and  Hoge,  1982) 
indicate  that  the  EM  bias  is  a  function  of  the  radar  wavelength.  The 
shorter  the  wavelength  the  smaller  the  value  of  EM  bias. 
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Standard  Deviation  of  EM  Bias  measurements  from  "best  fit" 
hyperspace  line*  All  of  the  63  possible  combinations  of  SWH, 
skewness,  kurtosis,  dominant  wavelength,  wind  speed,  and 
significant  slope  were  used  to  fit,  in  a  least  squares  sense, 
the  measurement  of  EM  Bias  to  a  hyperspace  line.  The  results 
are  displayed  in  the  order  of  highest  to  lowest  standard 
deviation.  The  dots  at  the  bottom  show  the  combination  used  to 
obtain  the  standard  deviation  plotted  above. 
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